
EVALUATION OF AN INFRARED 
HEATING SIMULATION OF 
A MACH 4.63 FLIGHT ON AN 
X-15 HORIZONTAL STABILIZER 

by Roger A. Fields and Andrew Vano 

Flight Research Center 
Fdwards, Calif. 


I%‘?0o2<c57f 

NASA TN D-5403 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • SEPTEMBER 1969 


I. Report No. 
TN D-5403 


2. Government Accession No. 


3. Recipient’s Catalog No. 


4. Title and Subtitle 

EVALUATION OF AN INFRARED HEATING SIMULATION 
OF A MACH 4. 63 FLIGHT ON AN X-15 HORIZONTAL 
STABILIZER 


5. Report Date 
September 1969 


6. Performing Organization Code 


7. Author(s) 

Roger A. Fields and Andrew Vano 


8. Performing Organization Report No. 
H-560 


9. Performing Organization Name and Address 
NASA Flight Research Center 
P. 0. Box 273 
Edwards, Calif. 93523 


10. Work Unit No. 

126-14-01-01-24 


1 I. Contract or Grant No. 


12. Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 
Washington, D. C. 20546 


13. Type of Report and Period Covered 
Technical Note 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

Temperatures recorded on the X-15 horizontal stabilizer during a Mach 4. 63 flight were simulated 
in the laboratory. A liquid-nitrogen evaporative cooler was used to cool the structure to a prelaunch 
condition; the heating was provided by an infrared heating system with closed-loop control. The simu- 
lated flight produced temperatures from approximately -50° F (228° K) to 750° F (672° K). The simu- 
lation was evaluated by comparing flight -measured temperatures with those measured during the 
simulation. 


IT. Key Words Suggested by Author(s) 
Flight heating simulation 


18. Distribution Statement 
Unclassified - Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

33 

$3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information, Springfield, 

Virginia 22151. 


J 


Page Intentionally Left Blank 


EVALUATION OF AN INFRARED HEATING SIMULATION OF A 

MACH 4. 63 FLIGHT ON AN X-15 HORIZONTAL STABILIZER 

By Roger A. Fields and Andrew Vano 
Flight Research Center 


SUMMARY 


Temperatures recorded on the X-15 horizontal stabilizer during a Mach 4. 63 flight 
were simulated in the laboratory. The simulated flight involved temperatures from 
approximately -50° F (228° K) to 750° F (672° K). A liquid-nitrogen evaporative cooler 
was used to cool the structure to a prelaunch condition; the heating was provided by an 
infrared heating system with closed -loop control. The heater was divided into areas 
called zones, with feedback from a control thermocouple in each zone. 

The simulation was evaluated by comparing simulation and flight temperature data. 
The thermocouples were grouped according to function or location on the horizontal - 
tail structure as follows: control, skin, web, and beam cap. The maximum average 
absolute temperature deviations over the complete simulated flight profile were 11 F° 

(6 K°) for the control thermocouples, 39 F° (22 K°) for the skin thermocouples, 30 F° 
(17 K°) for the web thermocouples, and 48 F° (27 K°) for the beam-cap thermocouples. 

It was found that accurate temperature simulation requires liberal and judicious 
use of shields between control zones and at major heat -sink boundaries, minimum size 
of heating control zones in which a mean zone temperature is programed for control, 
and careful heater design. 


INTRODUCTION 


The loads carried by the structure of a high-speed vehicle can be classified as 
(1) loads caused by aerodynamic forces and (2) loads caused by thermal stress in the 
structure itself (ref. 1). When strain gages are used to measure these loads, some 
method must be utilized to establish the response of the strain gage to its thermo - 
structural environment. 

One method of determining the thermostructural load effect of gage response to 
the flight environment is to observe the strain-gage output while duplicating the flight 
structure temperature distribution without aerodynamic loads. These responses may 
then be subtracted from the flight measurements to obtain aerodynamic flight load 
measurements. 

This report evaluates a method of simulating the structure temperature distribution 
resulting from aerodynamic heating. Only thermocouple data were used; no attempt 


was made to establish strain -gage responses. The method used infrared heating lamps 
with appropriate control equipment to match flight -measured temperature time histories 
on an X-15 horizontal stabilizer. The temperature range was approximately -50° F 
(228° K) to 750° F (672° K). The evaluation is based on a comparison of flight and 
laboratory test temperatures at selected points on the structure. 

The units used for physical quantities in this report are given both in U. S. Cus- 
tomary Units and the International System of Units (SI). Factors relating the two sys- 
tems are presented in reference 2; those used in this paper are presented in appen- 
dix A. 


TEST ARTICLE 


Description 

The X-15 airplane, a high-speed, high -altitude research vehicle, was used by the 
NASA Flight Research Center to investigate the many aspects of hypersonic flight. The 
aircraft specifications are presented in the tables of reference 3. 

The X-15 horizontal stabilizer served as a pitch and roll control and is located as 
shown in figure 1. The stabilizer is of conventional semimonocoque construction 



Figure 1.— X-15 airplane. 
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(fig. 2). Although Inconel-X is the primary material, the aft ribs and trailing -edge 
beam are a titanium alloy, and the leading -edge beam is stainless steel. The total 
exposed surface area is 51.7 square feet (4. 80 square meters). Additional structural 
details are presented in appendix B. 
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Figure 2.- X-15 horizontal-stabilizer structure. 


Instrumentation 


The horizontal stabilizer is instrumented internally with 121 thermocouples, and 
11 thermocouples are located on the outer surface of the torque box, as shown in fig- 
ure 3. Identification numbers are given in the figure only for the thermocouples spe- 
cifically discussed in this report. Thirty -gage (American Wire Gauge) chromel-alumel 
wire was used for all thermocouples. The instrumentation installation on the aft portion 

~TT ':u I x x-» Y 

Section A-A Section B-B Section C-C Section D-D Section E-E 


Section F-F Section G-G Section H-H Section I -I Section J-J 



O Upper and lower surfaces 
□ Surface centerline 
Beam or rib section 


Figure 3. - Thermocouple locations on X-15 horizontal stabilizer. Sections drawn double-size; number in 
parentheses denotes thermocouple number on lower surface. 
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of the stabilizer is shown in the photograph of figure 4. Strain gages were installed 
adjacent to the thermocouples on beams and ribs and at all thermocouples on the skin 
at section A -A of figure 2. Although strain-gage data are not included in this report, 
the gages were operated during simulation tests since they produce a small amount of 
local heating. 



E- 160 76 

Figure 4. - Instrumentation installation on X-15 horizontal stabilizer. (Aft panel skin removed.) 


FLIGHT DESCRIPTION 


The instrumented horizontal stabilizer was flown on the number 3 X-15 airplane 
during the flight from which the data presented in this report were obtained. The X-15 
was launched from the B-52 carrier aircraft at an altitude of 46,700 feet (14,230 me- 
ters) and a Mach number of 0. 82. Rocket-engine ignition occurred immediately after 
launch, and engine shutdown occurred about 90 seconds after launch. Immediately after 
shutdown, a horizontal -stabilizer load maneuver was initiated at a Mach number of 
approximately 4.5, an altitude of 81,000 feet (24,690 meters), and a dynamic pres- 
sure of 840 lb/ft 2 (402 hN/m 2 ). The maneuver consisted of a 90° left bank, a rapid 
increase in angle of attack to 11° and 3g normal acceleration, then a push-over to Og. 

Three other maneuvers during the flight were similar to the stabilizer load maneu- 
ver in that they were pullups to approximately 12° angle of attack and were performed 
at Mach numbers of 4. 0, 3.7, and 2. 5, dynamic pressures of approximately 635 lb/ft 2 
(304 hN/m 2 ), 615 lb/ft 2 (295 hN/m 2 ), and 645 lb/ft 2 (309 hN/m 2 ), and flight times of 
130 seconds, 150 seconds, and 220 seconds after launch, respectively. 
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The maximum Mach number, altitude, and dynamic pressure recorded during the 
flight were 4. 63, 84,400 feet (25,730 meters), and 1399 lb/ft z (670 hN/m z ), respec- 
tively. 


LABORATORY SIMULATION EQUIPMENT 


Heating System 

Infrared lamps were used to provide the radiant heat flux for the aerodynamic 
heating simulation. The lamps were mounted on polished stainless -steel reflectors of 
the same shape and contour as the stabilizer; the reflectors were then positioned 
5.5 ±0.5 inches (13. 97 ±1. 27 centimeters) above and below the stabilizer to form the 
primary heater. The top reflector assembly is shown in figure 5. The quartz heating 



Figure 5.— Top reflector assembly for heating horizontal stabilizer. 
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lamps were distributed on the reflector with the highest lamp densities over the main 
beam, leading -edge assembly, and trailing-edge beam. No attempt was made to pro- 
vide heating simulation around the stabilizer support, since this was deemed to be 
beyond the scope of this investigation. 

The infrared heater is controlled by a closed-loop system, as illustrated by the 
block diagram of figure 6. A control thermocouple at a particular specimen location, 
chosen to be representative of the zone, generates a millivoltage which is compared 
with the function-generator output; this output represents the simulation temperature 
time history programed for the control thermocouple location. If there is a differ- 
ence (error) between these two signals, the controller commands the ignitron power 
regulator to supply more or less power to the heater, maintaining the error near 
zero. 



Figure 6.- Block diagram of laboratory heating system. 


Typical power controller racks are shown in figure 7. Each rack contains two 
function generators which can program any of the six temperature controllers in the 
same rack. The regulator voltage for each channel and a local control board are also 
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E- 14609 

Figure Z— Laboratory heating system power controller rack. 
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displayed. The test is conducted from the central control console shown in figure 8. 
A selected number of temperatures are monitored on the strip-chart recorders and 
bar-graph display during the test. 



Figure 8.- Laboratory test control room. 


Additional control features of the heating system that improve the temperature 
control are rate, reset, and gain. The rate control modifies the command signal in a 
way that prevents temperature overshoot caused by the thermal capacity of the heating 
lamps. The reset control eliminates the temperature difference between command and 
test specimen caused by thermal losses. The gain control adjusts the width of the 
proportional control band. The setting for any one of these controls depends on the 
setting of the others and must be experimentally established prior to the data run for 
minimum control error. 


The temperature controller specifications are as follows: 


Null measuring accuracy ±0. 25 percent of full scale or 

15 microvolts, whichever 
is greater 

Repeatability 0.1 percent 

Response time 20 milliseconds 

Proportional -band range 0.3 percent to 30 percent for 

a 20 -microvolt range 


The specifications for the function generator are as follows: 


Time -base accuracy . 

Dead band 

Repeatability 

Potentiometer linearity 


1. 0 percent of elapsed time 
0.01 percent of full scale 
0.05 percent of full scale 
0 . 2 percent 
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The accuracy with which the flight temperature time histories were programed 
for the temperature controller was evaluated by comparing the output of the function 
generators with the flight data. This comparison established the programing to be 
within ±5 F° (±3 K°) of the flight -measured temperatures. 

X-15 launch was chosen to be zero heating test time. It was found that errors in 
synchronizing the function generators to this time caused temperature errors of ±5 F° 
to ±15 F° (±3 K° to ±8 K°). 


Cooling System 


The simulation at a particular data event must include all significant environmental 
conditions beginning from some data reference. For loads measurements in which 
strain gages are used, the simulation must be started at the last available zero-load 
condition; for the X-15, this was just prior to B-52 takeoff. Thus, the horizontal- 
stabilizer simulation had to start at ambient conditions, cool down and cold soak to the 
launch condition, and then proceed to the X-15 flight heating profile. 


A cooling system was constructed to provide the cool-down and cold-soak portion 
of the simulation. During a particular test, the cooler mixes ambient air with liquid 
nitrogen (LN 2 ), which evaporates and extracts the heat of vaporization from the ambient 

air, cooling it to the desired temperature. This air is then directed to the test speci- 
men by a system of ducts. 


The construction of the cooler is illustrated in figure 9; the side plate is removed 



Figure 9. - Liquid-nitrogen evaporative cooler for stabilizer tests. 


( Side plate removed.) 
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to expose the interior. A squirrel-cage fan drives ambient air through an expansion 
nozzle into an evaporation chamber. Bypass baffles force the flow in the chamber near 
the walls where it mixes with LN£ droplets sprayed from the LN 2 nozzles. The LN 2 

was supplied by a portable 150 -gallon (0. 568-meter^) cryogenic tank; the fluid was 
transferred by pressurizing the tank with a controlled LN 2 boiloff. 

Temperature was controlled for a given airflow by manually regulating LN 2 line 

pressure such that a particular downstream duct temperature followed a predetermined 
profile. 


DATA -RECORDING EQUIPMENT 


X-15 Pulse Code Modulation System 

Flight temperatures were measured by an onboard pulse code modulation (PCM) 
encoder; the encoded measurements were transmitted to a ground telemetry station 
during the flight, decommutated, and recorded on tape. The data were then reduced 
to temperatures by using a digital computer programed to convert system counts to 
degrees Fahrenheit on the basis of the chromel-alumel thermocouple table in refer- 
ence 4. The accuracy of the system has been demonstrated to be better than ±0. 3 per- 
cent of full scale, or ±7 F° (±4 K°) at 1500° F (1090° K) and ±10 F° (±6 K°) at -160° F 
(167° K). 


Loads Calibration Laboratory Data-Acquisition System 

The Loads Calibration Laboratory data -acquisition system consists of portable 
acquisition units located near the sensors and a central control system. Each port- 
able unit provides signal conditioning, multiplexing, and 12 -bit analog-to -digital 
conversion. The digital signal is relayed to the central control system and a high- 
speed digital computer shown in figure 10. The computer provides all necessary con- 
trol for the data acquisition and formats the incoming digital data from the acquisition 
units for recording on a 9-track magnetic tape. Voltages from up to 800 transducers 
(400 thermocouples, 320 strain gages, and 80 potentiometers) can be recorded in 
ranges varying from ±5. 0 millivolts full-scale with a resolution of 2. 5 microvolts 
to ±4. 0 volts with a resolution of 2. 0 millivolts. The transducer signals are com- 
mutated at a nominal sampling rate of 10 samples per second per channel. The data- 
acquisition system can record the input voltages to an accuracy of 0. 5 percent of 
range maximum and 0. 2 percent of range typical (±1. 6 F° (±0. 89 K°) at 700° F (645° K) 
and ±2.0 F° (±1. 11 K°) at -50° F (228° K) on 20-millivolt range). The largest contri- 
bution to temperature -measurement errors is made by the thermocouple itself. How- 
ever, since the same thermocouples are used in flight and ground simulation, this 
error in the simulation tends to cancel. The thermocouples are considered to be 
accurate to the larger of ±5 F° (±3 K°) or ±0. 75 percent of the temperature. 
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Figure 10.- Laboratory data-acquisition system central control. 
DEVELOPMENT OF HEATER 


Heating-Reflector Configuration 

The initial distribution of the lamps into zones (areas to be controlled by one 
thermocouple) and within each zone is shown in figure 11; all the quartz lamps in any 



Zone boundary 


Figure 11.- Initial lamp zoning for X-15 horizontal-stabilizer heating simulation. 
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A condition caused by a control thermocouple being at a higher temperature than programed because of heating from 
outside the subject zone. This causes the power to the zone lamps to remain off until the measured temperature becomes 
less than that programed. 


one zone are wired in parallel and controlled by the same regulator and therefore have 
the same voltage. The distribution of lamps within each zone helped provide the proper 
distribution of heat flux over the stabilizer surface. This configuration consisted of 
222 12 -inch 1000-watt lamps per reflector. The size of the control zones ranged from 
approximately 1. 0 square foot to 3. 9 square feet (0. 10 square meter to 0. 36 square 
meter), and the number of lamps per zone ranged from 11 to 26. 

The preliminary heating -test results indicated that more lamps were necessary to 
heat the leading -edge structure. In addition, the power required to heat the leading- 
edge zones overheated the adjacent zones at the control -thermocouple locations, causing 
a lack of adequate control in these adjacent areas. This problem, known as cross- 
talk, 1 was corrected by installing a leading -edge heater and an asbestos shield between 
the leading -edge zones and adjacent zones (fig. 12). 


E- 19433 

Figure 12.- Heater assembly for heating simulation on X-15 horizontal stabilizer. 
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Further heating tests indicated that it would be necessary to increase the number 
of control zones over the stabilizer surface. Heaters were also constructed for the 
root rib to simulate flight heating in the gap between the stabilizer and the fuselage; 
the root-rib heaters were controlled by internal thermocouples on the root -rib 
centerline. Figure 13 shows the root -rib -heater construction. A row of lamps was 
installed on the root -rib reflector above the stabilizer skin surface to reduce the 
primary heater end effects; these lamps were wired to be part of adjacent skin control 
zones. Ten of the 25 lamps located parallel to the root rib are 18 inches (46 centi- 
meters) long and are rated 1600 watts at 240 volts; the remainder are 12 -inch (30 -centi- 
meter) 1000 -watt lamps. 



Figure 13.— Horizontal-stabilizer root-rib heater. 
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The final lamp zoning configuration, composed of 36 temperature controllers and 
power regulators, is shown in figure 14. This configuration has 496 quartz lamps, 
and all were operated in the 0 to 480 -volt range. The leading -edge reflector lamps, 
indicated by the first row of lamps in the three leading -edge zones, are controlled by 
the zones on the lower reflector. Additional tests indicated that it was necessary to 
replace the 1000-watt lamps over the stabilizer main beam with 2000-watt lamps in 
order to obtain increased heating in this area. 
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Figure 14.- Final lamp zoning for X-15 horizontal-stabilizer heating simulation. 

The overall test setup in its final configuration is shown in figure 15. The power 
regulators are in the background; a movable data-acquisition site is shown to the 
right. 



Figure 15.- Overall heating simulation test setup. 
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Surface Heat Flux 


In-flight aerodynamic heating is a function of the difference between the recovery 
temperature and the surface temperature. Laboratory radiant heating is a function of 
the fourth powers of the lamp temperature and the surface temperature. Because of 
this fundamental difference in heating, it is important that lamps be distributed to ac- 
commodate this variation within a control zone, particularly in locations where there 
is a significant heat sink. A computer program was written during the test program 
to evaluate the surface flux distribution for a specified lamp configuration; the pro- 
gram is described in appendix C. 

Lamp configuration data from the horizontal -stabilizer simulation heater at 
section A -A of figure 2 were entered into the lamp flux program for analysis. Fig- 
ure 16 compares the net radiant flux required to maintain the flight temperature time 
history with that provided by the lamps as computed by the program for a time 75 sec- 
onds into the flight. The flight data points are calculated from the time -rate of change 
(slope) of flight temperature data and converted to a required net radiant heat flux 
necessary to match this slope by multiplying the specimen-absorbed heat flux by a 
correction factor to correct for surface reflectivity. 
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Figure 16.- Normal surface flux distribution on cross section of horizontal stabilizer. 


The lamp flux program calculated the flux distribution by using a lamp power level 
determined by matching the total power under the calculated curve to a curve faired 
through the flight data points for each zone to a tolerance of ±5 percent. The leading- 
edge zone was not considered, since the instrumentation in this area was inadequate. 

A comparison of the calculated curve and the flight data indicates that the distribution 
is well matched except in areas where there is substructure. In these locations, the 
structural discontinuity is not compensated for by a corresponding flux discontinuity; 
this flux discontinuity could be accomplished by placing an asbestos shield between 
lamp zones similar to that at the leading -edge zone. 
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End Effects 


Lamp end effects were the major source of temperature errors at the extremes of 
the specimen surface. The difficulty is demonstrated by curve 1 in figure 17, in which 
computer calculated flux distributions are plotted against distance from the specimen 
edge. The calculations are generated on a reflector/lamp configuration similar to that 
of the stabilizer. The 75-watt/inch (30-watt/centimeter) lamps were spaced 4.25 inches 
(10. 80 centimeters) above the specimen, 1.3 inches (3.3 centimeters) below a re- 
flector of reflectivity 0. 70, and at 1-inch (2. 54 -centimeter) intervals over a 100 -inch 
(254 -centimeter) length. Curve 1 represents a condition in which no attempt was made 
to correct the end effect, and the flux dropped to 56 percent of the 8. 1 Btu/ft 2 -sec 
(919. 3 hW/m 2 ) inboard mean value. The first attempt to improve this condition is 
shown as curve 2; the lamps were extended 10 inches (25.4 centimeters) beyond the 
surface boundary, which increased the terminal flux to 90 percent of the inboard mean. 


Distance from end of specimen, cm 



Figure 17.- Infrared heater end effects. 


By locating an end reflector of 0. 70 reflectivity at the surface boundary, the dis- 
tribution was further improved and the flux at the terminal point was increased to 
94 percent as shown in curve 3. Increasing the lamp density at the boundary by adding 
four lamps at 0.5 inch (1.27 centimeters), 1.5 inches (3.81 centimeters), 2.5 inches 
(6. 35 centimeters), and 3. 5 inches (8. 89 centimeters) did not relieve the end effect as 
presented by curve 4. If a single lamp at 0. 5 inch (1. 27 centimeters) is added to the 
end -reflector configuration of curve 3, the terminal flux is greater than the inboard 
mean and the problem is reversed (curve 5). 

The method chosen to reduce the end effect necessarily depends on the particular 
heating problem and the extent to which it must be reduced; the data presented indicate 
that an end reflector can be a useful tool in reducing this type of error. 
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Additional Problem Areas 


Precautions must be taken during the design of the reflector to insure that the re- 
flector shape does not change significantly enough during heating that the heating dis- 
tribution is altered or the lamps and specimen contact, or both. The thermal deforma- 
tions are caused by internal thermal stress, which causes reflector panel instability 
(buckling) and thermal expansion. The primary reflector system on the horizontal 
stabilizer was not designed to relieve thermal deformations, since it is relatively 
small and not subject to large temperature variations. A small amount of reflector 
buckling was observed during the test, but the deformations caused by these buckles 
were small enough in magnitude and mode to be considered negligible. 

As first designed, the leading-edge heater consisted of a stainless steel U-shaped 
reflector with lamps backed by a steel Unistrut beam supported at its extremes (i. e. , 
at the root and at the tip of the stabilizer). During heating, the thermal gradient 
through this heater caused the heater to deform toward the leading edge until it made 
contact. The design was modified by removing the Unistrut and cutting the heater into 
three sections, each with its own supports. 

It was necessary to insure that reflector/specimen alinement was maintained be- 
tween tests when the reflectors were removed for modifications. 


Cooling -System Changes 

Several cooling -system alterations were required during the preliminary test pro- 
gram. The dominant deficiency of the cooler was its inability to provide uniform 
temperature flow across the stabilizer span. This difficulty was caused by a non- 
uniform mixing of air and LN 2 in the evaporation chamber, which resulted from two 
factors: (1) a conical spray used to fill a rectangular section, and (2) flow reversal. 
The first difficulty was overcome by inserting rectangular plates along the nozzle 
centerline, thus causing the LN 2 to evaporate away from the nozzles. The evaporation 
chamber flow was improved by contouring the outside walls for a gradually increasing 
cross-section change. 


Test Configuration and Procedure 

The final heating configuration consisted of 36 control zones. The largest zone 
was approximately 2 square feet (0. 19 square meter). 

The flight -test condition of primary concern in the heating simulation was that from 
X-15 launch to approximately 200 seconds of free flight. The laboratory test procedure 
used to obtain this flight -test condition consisted of (1) cool-down and cold-soak, (2) 
transfer to heating mode, and (3) conduct the heating test. During the initial tests, an 
attempt was made to control the temperature of the ambient air that was directed over 
the simulation surface to a temperature profile taken from a standard atmosphere table 
(ref. 5) and the B-52 flight profile. During the cool-down, the heating system was 
programed to maintain the control thermocouple temperatures at the launch temperatures. 
The lamps did not operate until the temperature cooled down to the programed tem- 
perature. 
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This method of cool-down did not simulate aerodynamic flight conditions (such as 
dynamic pressure and Mach number) and so resulted in an unsatisfactory simulation. 
The procedure was changed to include a rapid cool-down to -50° F (228° K). As the 
internal structure temperatures approached launch temperatures, the air temperature 
was decreased to -70° F (217° K), which is the standard ambient air temperature at the 
launch altitude. The cooling profile is shown in figure 18 as duct air temperature data; 
duct 1 is farthest inboard. The cooling air was turned off as dictated by internal tem- 
perature, and the heating test was started. To prevent a sudden surge of power to the 
heaters, the function generators were reprogramed at the start, or zero time, to 
-100° F (200° K) and linearly increased to the flight temperature profile in about 5 sec- 
onds. 



Figure 18.- Cold-soak temperature time history. 
RESULTS AND DISCUSSION 


Some method of processing the large volume of data obtained in this evaluation be- 
came necessary in order to compare the flight -measured temperatures with those sim- 
ulated during the ground heating tests. The method chosen included the following steps: 
(1) The data were thinned to one sample per second; (2) flight and test simulation time 
histories were plotted to identical scales so that overlays could be constructed; (3) the 
temperature errors were measured from the overlays and tabulated for all channels at 
15-second intervals, and (4) the average absolute errors were established. 

The data were grouped as control, skin (other than control), beam-cap, and web 
thermocouples to provide a basis for comparison. Thermocouple 5, located on the 
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root rib as shown in figure 3, was included with the skin temperature data rather than 
with the control data because the lamps in that control zone did not operate. A com- 
bination of cross-talk and conduction heated the area sufficiently to maintain the tem- 
perature at the control location higher than that programed on the function generator. 
Typical plots of flight and simulation test temperatures for the four thermocouple 
groups as a function of time are shown in figures 19 to 22. 
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Figure 19.- Typical plots of flight and simulation test control thermocouple temperatures versus time. 
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(a) Thermocouples 1,29,37, 54. 


Figure 20.— Typical plots of flight and simulation test skin temperatures versus time. 
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(b) Thermocouples 98, 99, 116, 117. 


Figure 20.- Concluded. 
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Figure 21.- Typical plots of flight and simulation test beam-cap temperatures versus time. 
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Figure 22.- Typical plots of flight and simulation test web temperatures versus time. 
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Figure 23 shows time histories of absolute temperature deviations for each of the 
four groups of thermocouples. Of all thermocouples included in each group, 25, 50, 
75, and 95 percent had deviations equal to or less than the respective curves shown in 
the figure. The control and skin thermocouple data temperature deviations increased 
during the portion of the test that corresponded to the times of maximum heating rate 
and maximum temperature. The beam -cap (heat sink) and web (internal) thermocouple 
data temperature deviations increased generally throughout the test as a cumulative 
effect of control and skin thermocouple errors. 


Percent of data equal to or 
less than this curve 


o 

25 

□ 

50 

A 

75 

o 

95 





-60 

“50 

-40 

-30 

-20 

-10 

.0 


Absolute 

temperature 

deviations, 

K° 


Figure 23.- Absolute deviations between flight and simulation test temperatures. 
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The averages of the absolute deviations between the simulation test and flight 
temperatures for each of the four thermocouple groups are shown plotted as a function 
of time in figure 24. The maximum absolute average temperature deviations for the 
control, skin, beam-cap, and web thermocouples are 11° F (6° K), 39° F (22° K), 48° F 
(27° K), and 30° F (17° K), respectively. The small temperature deviations of the con- 
trol thermocouples are caused by inaccuracies in heating control equipment and errors 
in the program timing of the function generators. 



Time, sec 


Figure 24.- Comparison of absolute average deviations between heating simulation and flight - 
measured temperatures for the four groups of thermocouples. 

In several areas of the stabilizer, temperature deviations were large. Most of 
these problem areas are typified by thermocouple 54 temperatures. The simulation 
temperature profile for thermocouple 54 (fig. 3) as plotted in figure 25 is similar to 
the temperature profile programed for the heating zone that included thermocouple 54. 
The flight temperature profile at this location is considerably different. This situation 
could have been improved by programing the mean measured or calculated temperature 
for the entire control zone. The variation of temperatures in a control zone depends 
on several factors, including the control -zone size and variations in the structure and 
flow conditions. These temperature variations can be reduced by increasing the 
number of control zones and decreasing the zone area. The number of control zones 
will be limited by the available control equipment (as in the test discussed herein), by 
the size of the lamps, and by the problems associated with control cross-talk. The 
need for shields at control -zone boundaries increases as the size of the zone decreases. 
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Figure 25. — Temperature variation in the heating zone that includes thermocouple 54. 

Generally, the temperatures on the main-beam caps were lower during the simu- 
lation than during flight (fig. 21). However, the web temperatures at these same 
locations were either higher or would be higher if the cap temperatures were correct 
(fig. 22). The discrepancy is attributed to the fact that the thin web is heated pri- 
marily by internal radiation rather than conduction. The skin immediately adjacent 
to the beam cap is overheated since it receives essentially the same flux as the cap; 
it therefore reradiates more to the web during the simulation than during the flight, 
causing the web to overheat. This condition further substantiates the need for shields 
at the boundaries of heat sinks. 

The leading -edge beam has a web three times as thick as the main beam and there- 
fore is not as subject to heating by internal radiation. The simulation temperatures 
on this web were found to be close to or slightly lower than the flight temperatures. 

Temperature deviations at other locations are an accumulation of errors due, in 
part, to the data-acquisition system, the thermocouples, the reflector design, struc- 
tural discontinuities, and restricting the programing and control of the temperature 
profile to a single location in a zone. 


CONCLUDING REMARKS AND RECOMMENDATIONS 


Temperatures recorded on the X-15 horizontal stabilizer during a flight to Mach 
4. 63 were simulated in the laboratory. The overall simulation was considered to be 
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good; maximum absolute average temperature deviations were 11 F° (6 K°) for the con- 
trol thermocouples, 39 F° (22 K°) for the skin thermocouples, 30 F° (17 K°) for the web 
thermocouples, and 48 F° (27 K°) for the beam-cap thermocouples. The temperatures 
at the control -thermocouple locations were close to those programed as recorded during 
flight. Errors at these locations were primarily inherent in the heating control equip- 
ment and the program timing of the function generators. 

Temperature errors at other locations resulted from a number of factors. The 
overall zone temperature simulation could have been improved if the control thermo- 
couples had been programed with a mean control zone temperature rather than the 
temperature time history at the location of the control thermocouples. Judicious 
location of shields between control zones and at significant structural discontinuities 
would help to eliminate cross-talk and improve local heat -flux distribution. Also, zone 
sizes should be kept to a minimum. 

Reflector design (lamp arrangement) is an important requirement for accurate 
temperature control. A systematic method to derive this design is necessary. The 
first step must establish the design criterion; a calculated flux distribution from the 
flight -measured data (or other data to be simulated) can serve as this criterion. At- 
tempts can then be made to obtain lamp distributions to match this flux distribution 
with appropriate discontinuities (i.e. , shields) at major heat -sink boundaries. A 
computer program which calculates the lamp flux distributions is a valuable tool for 
this purpose. Heater end effects must also be considered as a part of this design. 


Flight Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif., May 15, 1969. 
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APPENDIX A 


CONVERSION OF U. S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Con- 
ference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 2). 
Conversion factors of the units used herein are given in the following table: 


Physical quantity 

U. S. Customary 
Unit 

Conversion 

factor* 

SI Unit 

Heat flux 

Length 

Pressure 

Temperature 

Volume 

Btu/(ft 2 -sec) 

ft 

in. 

lb/ft 2 

°R = °F + 460 
gal 

1. 1349 x 10 4 
0.3048 
2. 54 
0.4788 
0.556 

3.785 X 10 " 3 

W/m 2 

m 

cm 

hN/m 2 

°K 

m 2 


*Multiply value given in U. S. Customary Unit by conversion 
factor to obtain equivalent values in SI Unit. 


Prefixes to indicate multiple of units are: 


Prefix 

Multiple 

centi (c) 
hecto (h) 

10 "2 
10 2 
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APPENDIX B 


HORIZONTAL -STABILIZER STRUCTURE 


The leading -edge beam, main beam, trailing-edge beam, and skins of the X-15 
horizontal stabilizer (fig. 2) form a two-cell torque box. The ribs are oriented per- 
pendicular to the main -beam centerline. The skin is made from a constant -thickness, 
0.050-inch (0. 127 -centimeter) Inconel -X sheet. The Inconel-X leading edge is designed 
to act as a heat sink and is constructed in five segments to relieve thermal stresses. 

The leading-edge beam is a U-channel made from stainless steel 0.093 inch (0.236 centi- 
meter) thick to station 37. 625 and tapered to 0. 050 inch (0. 127 centimeter) at the tip. 

The main beam is constructed entirely of Inconel-X. The corrugated web, 

0. 032 inch (0.081 centimeter) thick, is double thickness from the root to station 
26.875. The main-beam caps taper in width from 6.0625 inches (15.399 centimeters) 
at the root to 1. 875 inches (4. 762 centimeters) at station 37. 625 and to 1. 625 inches 

(4. 128 centimeters) at station 53. 750. The remainder of the outboard cap is a constant 

1. 625 inches (4. 128 centimeters) wide. The cap thickness along the centerline of the 
beam tapers from 0. 300 inch (0. 762 centimeter) at the root to 0. 096 inch (0. 244 centi- 
meter) at station 53. 750 and remains constant outboard. The beam height tapers from 
3. 976 inches (10. 099 centimeters) at the root to 1. 223 inches (3. 106 centimeters) at the 
tip. The main-beam caps extend inboard of the root rib to form the top and bottom of 
the structural box (torque box) that provides support for the stabilizer. 

The trailing-edge beam is made from a titanium alloy. It is a solid beam 
0.448 inch (1. 138 centimeters) wide that tapers in height at the aft edge from 0. 768 inch 
(1. 951 centimeters) inboard to 0. 178 inch (0. 452 centimeter) outboard. 

The ribs forward of the main beam are 0.040 inch- (0. 102 centimeter-) thick 
Inconel-X; those aft of the main beam are 0. 050 inch (0. 127 centimeter) titanium alloy. 
The ribs are spaced 5.375 inches (13. 652 centimeters) apart. 
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APPENDIX C 


LAMP FLUX COMPUTER PROGRAM 


The FORTRAN IV subset E computer program described in this appendix was 
written to calculate the normal surface heat flux distribution (Btu/ft 2 -sec) from a 
parallel nonuniform radiant heater system. The program is intended to serve as a 
method of predicting the effectiveness of a particular infrared heating system design 


Heater Model 

The heater -specimen configuration shown in the sketch below was chosen to idealize 
the typical test setup. N -lamps are positioned parallel to and D distance from the flat 
specimen surface; the lamp spacing is defined by the array Y(I), I = 1,N which is 
referenced to X = 0. The lamps are considered to be infinite in length (dimension per- 
pendicular to the plane of the figure) and are backed by a reflector of reflectivity U1 


Y(N) 

Y(N-l) 

YIN -2) 

YIN -3) 

Yll) 

YM) 

YI3) 

YI2) 

Yll) 


Shield U2 



(defined as the ratio of reflected to incident radiation), which is parallel to the specimen 
surface and infinite in both dimensions . A shield of reflectivity U2 is positioned at 
the reference plane X = 0. The variables are further defined in the following table: 


Variable 

D 

R 

U1 

U2 

M 

N 

X0 

XINC 

Y(I) 

WATT 


Definition 

Lamp to surface dimension, real 
Lamp to reflector U1 dimension, real 
Reflector reflectivity, real 
Shield reflectivity, real 

Number of X 


calculation locations, integer, M = 


+ 1 


Number of lamps, integer 
X value at start of calculations, real 
Calculation interval, real 
Shield to i th lamp dimension (array), real 
Lamp output: watts/inch, real 
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The calculation range starts at X = XO and moves to XMAX in M number of XINC 
increments as defined by the equation. The lamp output is defined by the parameter 
WATT and represents the lamp rating in watts per inch at the selected power level. 

The idealization suffers from a parallel lamp-reflector -surface restriction and a 
uniform longitudinal lamp assumption; however, by choosing suitable geometry and 
reflector/shield configurations, a large spectrum of problems can be solved. 


Method of Calculation 

A superposition method is used to calculate the specimen surface flux. The flux at 
a particular X is calculated by using an inverse flux/distance -to -lamp relationship 
which is corrected for angle of incidence. The flux is then summed for each lamp by 
using the program listing below and the equation on the following page. 


Program listing 

OIMFNSION Y< 500) 

1 RFAD(i,2)D*R*Ul,U2,M,N,XO,XINC,WATT t (YU ) *1 = 1 tM 

2 FORMAT ( 2F6.2, ?F5.2,2I4,?F6.2,FR.?/{ 12F6. I ) ) 

IF{0-999.<39) 3,4 ,4 

3 R 6 AO ( 1,5) 

WR T TF ( 3,5) 

5 FORMAT ! 50H 
WR TTF( 3,12) 

12 FORMAT! 54H0- CALCULATION PARAMETERS 0*R»U1,U?,M,N,XG,XINC,WATT,YS) 
WPITEt ?,2)D,R»U1,U2,M,N,X0,XINC,WATT , ( Y ( I ) , I = 1 , N ) 

WR ITF! 3,11) 

11 FORMAT ( 5 7H0 X(rN) RTU/F T2-SEC DIP REF S HD RFSH 

10) 

PEAL K0*K1,K?,K3,K4,K5 

K0=D**2 

K 1=D + 2*R 

K2=U1*K1 

K3=U?*0 

K4=U?*K2 

K5=K 1**2 

6 00 7 1=1, M 
X=X0 + ( I- 1 ) *X I NC 

sumw~o 

SUM W 1= 0 
SUMW2=C 
SUMW7=0 
SUMW4=C 
8 no 9 J = 1 , N 
V2=Y(J)+X 
v 3 = y ( j )-y 
V4= V 2* * 2 
V*=V 3**2 

Wl=0/( 2*3. 1416*! K0+V5) ) 

W2=K2/! 2*3. 1416*(K5+V5) ) 

W3=K3/t 2*3. 1415*! KC+V4) ) 

W4=K4/! 2*3. 1 4 1 6* ! K 5+ V4 ) ) 

W=W1+W2+W3+W4 
Sl'MW=9UMW-»-W 
SUMW1=SUMW1+W1 
SUMW?=SUMW2+W2 
SUMW3= SUMW3 + W3 
q SUMW4=^UMW4+W4 
DIR=SUMW1 /SUMW 
RFF=SUMW2/SUMW 
SH0=SUMW3/SUMW 
RFS0=SUMW4/SUMW 
PTU=0 . 1365*SUMW*WATT 

7 WR ITE! 3, 1C) X,BTU,DIR ,REF ,SH0,RESQ 

10 FORMAT!F8.2,F11.3, • »,2 P4F9.2) 

GO TO 1 

4 WRITE! 3, 13) 

13 FORMAT ! ' CEND CALCULATIONS') 

STOP 

FNO 
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FLUX 


( Btu 
l ft 2 -sec 


) -(NF) ™ ATT * sh— * 

'] v ' jti Id 2 + (Y(j 


(Y(j) - X)2 


U1(D+ 2R) 


U2(D) 


U1(U2)(D + 2R) 


(D+ 2R) 2 + (Y(J) - X) 2 D 2 +(Y(J)+X) 2 (D+ 2R) 2 + (Y(J) + X) 2 . 


The first term of the equation represents the normal surface flux due to direct 
radiation; the second term, that due to radiation reflected from the parallel reflector 
Ul. The third term represents the normal surface flux caused by radiation from the 
vertical X = 0 reflector U2, and the last term calculates the radiation coming from 
both reflectors Ul and U2. The program is written to calculate flux in U. S. Custom - 
ary Units; the constant (0. 1365) must be changed if flux in other units is desired. 


Input 


To use the program a punched data deck must be generated according to the fol- 
lowing format: 


Parameter 

D 

R 

Ul 

U2 

M 

N 

XO 

XINC 

WATT 

Columns 

1-6 

7-12 

13-17 

18-22 

23-26 

27-30 

31-37 

37-42 

43-50 

Decimal 

4 

10 

15 

20 

INT 

INT 

34 

40 

48 


Card 1: Control information 


Parameter 

Columns 

Decimal 


Y(l) 

Y(2) 

Y(3) 

Y(4) 

Y(5) 

Y(6) 

Y(7) 

Y(8) 

Y(9) 

Y(10) 

Y(ll) 

Y(12) 

1-6 

7-12 

13-18 

19-24 

25-30 

31-36 

37-42 

43-48 

49-54 

55-60 

61-66 

67-72 

5 

11 

17 

23 

29 

35 

41 

47 

53 

59 

65 

71 


Lamp array cards (typical) 


Parameter 

Columns 


Any title 

2-50 

Title card 


Enter 

Columns 

Decimal 

Enter 

Columns 

Decimal 


999.99 

0.0 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

1-6 

7-12 

13-17 

18-22 

23-26 

27-30 

31-36 

37-42 

43-50 

4 

10 

15 

20 

INT 

INT 

34 

40 

48 


0.0 


1-6 

4 

Stop cards 


This data deck is then used with the lamp flux object deck. As many problems (data 
sets) as required may be processed at one time by stacking groups of the control in- 
formation, lamp array, and title cards in order. The lamp array is limited to 500 lamps 
per problem. The two-card stop is required to terminate the calculation. 
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